Mean shear flows generated by nonlinear resonant Alfven waves 
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In the context of resonant absorption, nonlinearity has two different manifestations. The first is the 
reduction in amplitude of perturbations around the resonant point (wave energy absorption). The 
second is the generation of mean shear flows outside the dissipative layer surrounding the resonant 
point. Ruderman et al. [Phys. Plasmas 4, 75 (1997)] studied both these effects at the slow resonance 
in isotropic plasmas. Clack et al. [Astron. Astrophys. 494, 317 (2009)] investigated nonlinearity 
at the Alfven resonance, however, they did not include the generation of mean shear flow. In this 
present paper, we investigate the mean shear flow, analytically, and study its properties. We find that 
the flow generated is parallel to the magnetic surfaces and has a characteristic velocity proportional 
to e 1 / 2 , where e is the dimensionless amplitude of perturbations far away from the resonance. This is, 
qualitatively, similar to the flow generated at the slow resonance. The jumps in the derivatives of the 
parallel and perpendicular components of mean shear flow across the dissipative layer are derived. 
We estimate the generated mean shear flow to be of the order of lOkms -1 in both the solar upper 
chromosphere and solar corona, however, this value strongly depends on the choice of boundary 
conditions. It is proposed that the generated mean shear flow can produce a Kelvin-Helmholtz 
instability at the dissipative layer which can create turbulent motions. This instability would be an 
additional effect, as a Kelvin-Helmholtz instability may already exist due to the velocity field of the 
resonant Alfven waves. This flow can also be superimposed onto existing large scale motions in the 
solar upper atmosphere. 

PACS numbers: 52.25.Fi; 52.30.Cv; 52.35.-g; 52.35.Bj; 52.35.Mw 



I. INTRODUCTION 

It has long been established, from observations, that 
the solar corona is highly structured and inhomogeneous 
with temperatures of the order of 10 6 K. The solar 
corona is filled with a large number of discrete magnetic 
loops (coronal arcades) and there is an abundance of ob- 
servational evidence showing that magnetohydrodynamic 
(MHD) waves propagate in, and across, these loops (see, 
e.g. Refs. In order to sustain such high tem- 

peratures whilst combatting optically thin radiation and 
thermal conduction there must exist some mechanism(s) 
acting as a source of steady heating. 

The last few decades saw a multitude of models pro- 
posed to tackle the complicated problem of coronal heat- 
ing, such as heating by waves (e.g. resonant absorption, 
phase mixing) or heating by magnetic relaxation (e.g. 
reconnection) . There is an increasing consensus that all 
these processes act simultaneously to different degrees of 
efficiency throughout the solar corona. The agreement 
comes as observational data confirms the abundance of 
MHD wave propagation in the solar corona, along with 
the complexity of the magnetic configuration enabling 
magnetic heating to occur (see, e.g. Refs. (5|-[Toj). 

Resonant absorption of Alfven waves in coronal loops 
was first suggested by Ionson(author?) [11] as a non- 
thermal heating mechanism of the corona and since then 
resonant absorption has been studied intensively by many 
authors in the linear regime for a variety of applications 



(see, e.g. Refs. [12|]-[25|). A new understanding of the 



process of resonant absorption became available after the 



paper by Ruderman et al. (author?) (2(| which was the 
first analytical study on the nonlinear aspect of resonant 
absorption. Not only was it shown that resonant ab- 
sorption of slow waves was an inherently nonlinear phe- 
nomenon, they also showed that a mean shear flow is 
generated outside the dissipative layer. By their calcu- 
lations, however, they still found the generated flow was 
much too large compared to the observed velocities. On 
the other hand, the authors did note that their results 
should be used with caution in the solar atmosphere as 
some of their assumptions were not fully realistic for that 
environment. 

From the coronal heating point of view, the slow reso- 
nance studied by Ruderman et al. (author?) [26[ was not 
expected to contribute significantly as the energy stored 
in slow waves is much less than the required energy to 
compensate the losses and it is difficult for slow waves to 
reach the corona, as they become shocked as they climb 
due to density stratification. From this point of view, 
Alfven waves and Alfven resonance are of much greater 
interest as estimations show that the energy carried by 
Alfven waves is much higher and they can reach the 
corona without significant damping. The validity of non- 
linear resonant Alfven waves (under coronal conditions) 
was studied in great detail by Clack et al. (author?) [27l ]. 
where they showed that in coronal plasmas the nonlin- 
ear addition to the result found in linear MHD is so 
small that the linear approach can be used with great 
accuracy. Essentially this investigation clarified the up- 
per limit in which linear theory is applicable to resonant 
Alfven waves. 

Studies have been carried out to investigate the prop- 
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erties of shear flows generated by velocity field of Alfven 
waves, however, nearly all of these have been numerical 
due to analytical complications when considering non- 
linearity, turbulence and resonant absorption simultane- 
ously. These studies have found that shear flows could 
give rise to a Kelvin-Helmholtz instability at the nar- 
row dissipative layer (see, e.g. Ref. 28]). This instabil- 
ity can drive turbulent motions and, in turn, locally en- 
hance transport coefficients which can alter the efficiency 
of heating (see, e.g. Refs. [28|-[30l]). None of these in- 
vestigations have studied the generation of mean shear 
flow at the Alfven resonance by nonlinear interactions. 
The generation of mean shear flow can supply additional 
shear enhancing turbulent motions. 

Recent advancements in the understanding of nonlin- 
ear Alfven resonance has inspired us to study of the gen- 
eration of mean shear flows at the Alfven resonance. In 
the present paper we will derive the equations describing 
the generated mean shear flow outside the Alfven dissipa- 
tive layer and estimate the magnitude of the shear flow. 
We already know, from Ofman and Davila(author?) 
[28| , that the plasma velocity at the dissipative layer may 
be reduced by the turbulent enhancement of the dissipa- 
tive parameters, implying that for a given heating rate, 
the wave amplitude is reduced compared to the linear 
case. This means that any result we produce must be re- 
duced when considered for the solar corona since turbu- 
lent motions are likely to be present and this will reduce 
the mean flow speed. 

The paper is organised as follows. In the next sec- 
tion, we introduce the governing equations and discuss 
the main assumptions. In Sect. Ill we recall some previ- 
ous results and derive the solution for the equation gov- 
erning the mean shear flow outside the dissipative layer. 
In Sect. IV we derive the governing equations for the 
mean shear flow inside the dissipative layer. Section V 
gives the jumps in the derivatives of the mean shear flow 
velocities across the dissipative layer and we estimate the 
magnitude of these jumps for conditions typical in the 
solar upper chromosphere and corona. In Sect. VI we 
draw our conclusions and discuss the results. Once the 
governing equations of the generated mean shear flow are 
found, the nonlinear theory of resonant Alfven waves is 
complete. 



II. GOVERNING EQUATIONS AND 
ASSUMPTIONS 

To mathematically study the mean shear flow gener- 
ated at resonance we use the visco-resistive MHD equa- 
tions. In spite of the presence of dissipation, we use the 
adiabatic equation as an approximation to the energy 
equation. Numerical studies by Poedts et al. (author?) 
[3l| in linear MHD have shown that dissipation due to 
viscosity and finite electrical conductivity in the energy 
equation does not, significantly, alter the behaviour of 
resonant MHD waves in the driven problem. 



When the product of the ion (electron) gyrofre- 
quency, lu^, and the ion (electron) collision time, Tj( e ), 
is much greater than unity (as in the solar corona) 
the viscosity and finite electrical conductivity become 
anisotropic. The parallel and perpendicular components 
of anisotropic finite electrical conductivity only differ by 
a factor of 2, therefore, we will consider only one of them 
without loss of generality. The anisotropic viscosity is 
given by Braginskii's viscosity tensor (see, e.g. Braginskii 
|32|). The components of the viscosity tensor that remove 
the Alfven singularity are the shear components, however 
these components (in their form) act as the isotropic vis- 
cosity. Indeed, Clack et aZ.(author?) (2?| showed that 
one can interchange isotropic viscosity and shear viscos- 
ity when studying the Alfven resonance. Another pos- 
sible additive to our model might be dispersive effects 
caused by Hall currents, however Clack et al. (author?) 
[27j showed that this effect does not alter the MHD waves 
in the vicinity of the Alfven resonance. 

The dynamics of waves in our model is described by 
the visco-resistive MHD equations 
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In Eqs. (JT|) — C|4j , v and B are the velocity and magnetic 
induction vectors, P the total pressure, p the plasma 
pressure, p the density, rj the coefficient of magnetic dif- 
fusivity, V = rji/p the coefficient of shear viscosity over 
density, V-Si the shear component of Braginskii's viscos- 
ity tensor, 7 the adiabatic exponent, and /io the magnetic 
permeability of free space. 

We adopt Cartesian coordinates x, y, z and limit our 
analysis to a static background equilibrium (vo = 0). 
We assume that all equilibrium quantities (terms with 
subscript '0') depend on x only. The equilibrium mag- 
netic field, Bo, is unidirectional and lies in the yz-planc. 
The equilibrium quantities must satisfy the condition of 
total pressure balance, po + B^/(2p ) — constant. 

For simplicity we assume that the perturbations of all 
quantities are independent of y (d/dy = 0). We note 
that since the magnetic field is not aligned with the z- 
axis, an Alfven resonance can still exists. Even though 
the Alfven resonance is governed by a linear equation, 
we must consider nonlinear effects to obtain the second 
manifestation of nonlinearity - mean shear flows. In the 
linear theory of driven waves all perturbed quantities os- 
cillate with the same frequency, u), which means that 
they can be Fourier-analysed and taken to be propor- 
tional to exp(i[fcz — Lot]). Solutions are sought in the 
form of propagating waves and all perturbations in these 
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solutions depend on the combination 6 — z — Vt, rather 
than z and t separately, with V = cu/k. In the context of 
resonant absorption the phase velocity, V, must match 
the projection of the Alfven velocity, va, onto the z-axis 
when x — x a where x a is the resonant position. To define 
the resonant position mathematically it is convenient to 
introduce the angle, a, between the z-axis and the direc- 
tion of the equilibrium magnetic field, so that the com- 
ponents of the equilibrium magnetic field can be written 
as B 0y — B sin a, B 0z = B cos a. The definition of the 
resonant position can now be defined mathematically as 



V = Va (% a ) cos a , 



(5) 



where va is the Alfven speed defined as va = 
Bof(noPo) 1 ^ 2 ■ In addition, we introduce the squares 
of the sound and cusp speeds as c| = 7P0/P0 an d 
c t = c s v a/(. c s + v a)> respectively. In what follows we 
can take x a = without loss of generality. The pertur- 
bations of the physical quantities are defined by 



p = p Q + p, p = p Q +p, 
B b b 2 
Mo 2^io ' 



B = B„ 



P 



P 



(6) 



where P is the perturbation of total pressure. 

The dominant dynamics of resonant Alfven waves, in 
linear MHD, resides in the components of the perturbed 
magnetic field and velocity that are perpendicular to the 
equilibrium magnetic field and to the x-direction. This 
dominant behaviour is created by an x^ 1 singularity in 
the spatial solution of these quantities at the Alfven reso- 
nance (see, e.g, Sakurai et al. [l7|, Goossens et al. [33|): 
these variables are known as large variables. All other 
variables are known as small variables. 

To make the mathematical analysis more concise, we 
define the components of velocity and magnetic field that 
are in the j/z-plane and are either parallel or perpendic- 
ular to the equilibrium magnetic field: 



w 



v — w 
by - b z 



sin a 
cos a 



cos a 
sin a 



(7) 



where v, w, b y and b z are the y- and z-components of the 
velocity and perturbation of magnetic field, respectively. 

Let us introduce the characteristic scale of inhomo- 
geneity, 1-^. The classical viscous Reynolds number, R e , 
and the magnetic Reynolds number, i? m , are defined as 



Re. = 



Vh 



inh 



Rn 



VL 



V 



(8) 



These two numbers determine the importance of viscosity 
and finite electrical conductivity. We introduce the total 
Reynolds number as 



1 

R 



1 

Re 



1 

Rn, 



(9) 



The aim of this paper is to study the generation of a 
mean shear flow outside the dissipative layer due to the 
nonlinear behaviour of driven resonant Alfven waves in 
the dissipative layer. We are not interested in the effects 
of MHD waves that have large amplitude everywhere and 
require a nonlinear description in the whole space. We 
focus on waves that have small dimensionless amplitude 
e <C 1 far away from the ideal Alfven resonant point. 

In nonlinear theory, when studying resonant behaviour 
in the dissipative layer, we must re-scale the dissipative 
coefficients (see, e-g^ Ruderman et al. [2(|, Clack et al. 

Ballai et al. [34| and Clack and Ballai (35|) so that 
dissipation is of the same order as nonlinearity, 



(10) 



It is easy to estimate the nonlinearity parameter, which 
is the ratio of the largest nonlinear and dissipative terms 
and is obtained to be ei? 2 / 3 (see Clack et al. [27]). This 
implies that if ei? 2 / 3 1 linear theory holds, whereas 
if ei? 2 / 3 > 1 then nonlinearity is important and linear 
theory breaks down. 

In linear theory all perturbed quantities are harmonic 
functions of 9, therefore their mean values over a period 
vanish. On the other hand, in nonlinear theory the per- 
turbed variables can have nonzero mean values as a result 
of nonlinear interaction of different harmonics. Let us in- 
troduce the mean value of a function f(6) over a period 
L as 



(/) = T / /(*) d0 - 



It directly follows from Eq. |T]) that 

(pu) = (b x ) = 0. 



(11) 



(12) 



We can always define the background state in such a way 
that the mean values of density, pressure and magnetic 
field vanish: 



(p) = (p> = (6|i) = (6x>=0. 



(13) 



This is not possible for the velocity, since we assume a 
static equilibrium. It is convenient, therefore, to divide 
Dm and v± into mean and oscillatory parts. Using the 
Reynolds decomposition we can write 



U\\ = {v\\), U± = {v±), 

v\\ = ^ 1 1 — 1 1 , v± = v± — U± 



(14) 



The quantities U\\ and U±_ describe the mean flow parallel 
to the magnetic surfaces. The mean flow is generated by 
the nonlinear interaction of the harmonics in the Fourier 
expansion of the perturbed quantities with respect to 6. 
We can rewrite Eqs. (HJ) — (J3J) m the scalar form as 
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b z db± 
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i?ou± cos a) = 



d(ub±) d(wb±) 



dx 

dv± , dvi 



de 
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w 
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dx 
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P dx ^ dx 



2B b 



(17) 



(18) 
(19) 



(20) 



^-"'W-^' (21) 



(22) 



(23) 



The equations for U\i and U± are obtained by averaging 
Eqs. (fT7|) and (JTSJ) , respectively, and then dividing by p, 
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(25) 



Here V = (d/dx,0,d/de), u is the x-component of the 
velocity, and w — v\\ cos a — v±_ sin a is the oscillatory 
part of w. 

Equations (fl5|) - ([23)l will be used in the following sec- 
tions in order to calculate the mean flow that is generated 
outside the dissipative layer by the resonant waves. 



III. OUTER SOLUTION 

To calculate the mean shear flow generated outside the 
dissipative layer we have to recall some results found by 
Clack et al. (author?) [27J. They used the method of 
simplified matched asym ptot ic expansions (see, e.g., Bal- 
lai et al. [34[ and Nayfeh [36j) in order to derive the equa- 
tion governing resonant Alfven waves in the dissipative 
layer. The method consists of finding the so-called outer 
and inner expansions and matching these expansions in 
overlap regions around the dissipative layer. This nomen- 
clature is well adopted for our scenario. The outer expan- 
sion corresponds to the solution outside the dissipative 
layer and the inner expansion refers to the solution inside 
the dissipative layer. The solution in the outer region is 
represented by asymptotic expansions of the form 



f = ef 



(1) _j_ £ 3/2j(2) 



(26) 



where e is the dimensionless amplitude of perturbations 
far from the dissipative layer and / represents any per- 
turbed quantity with the exception of vn and v±. We 
shall show that resonant Alfven waves create a shear flow 
with an amplitude proportional to e 1 / 2 outside the dissi- 
pative layer. As a consequence, we expand Vh and v±_ in 
the form of Eq. (f2"6"| , and U\\ and U± in the form 



£ V2 [7 (0) +eC/ (l) + e 3/2 [/ (2) + 



C/ ± = £ 1 /2 t/ (0) +eC/ (l) +e 3/2 [/ (2) + 



(27) 



Clack et al. (author?) 27] found that the substitution 
of Eq. (26]) into Eqs. ([T5]) - d23j) leads, in the first or- 
der approximation, to a system of linear equations for 
the variables with the superscript '1'. All variables can 
be eliminated in favour of vS 1 ' and P^\ leading to the 
system 



V- 



de 



F 



dx 



V- 



dx 



PoD A 



de ' 



where 



F 



p D A D c 



V 4 - V 2 {v\ + c|) + v\ c| cos 2 a ' 



(28) 



(29) 



5 



Da = V 2 — v\ cos 2 a, 



Dc = K + c s) - c 



2 cos 2 a) 



(30) 



The quantities Da and Dc vanish at the Alfven and 
slow resonant positions, respectively. As a result these 
two positions are regular singular points for the system 
(|28[) . The remaining variables can be expressed in terms 
of «W an d p(i) as , 



__.(!) = Vsma p(1) 



PoDj 



Ucg cos a 



P (1) , (31) 



/-n BoCosq /-n (i) Bocosasina 



P (1) , (32) 



where VJ^ and V± are constants and the superscripts '-' 
and '+' refer to x < and x > 0, respectively. 

Continuing in the same way we can find the solutions 
of the subsequent higher order approximations. At each 
step in the scheme of approximations we obtain equations 
with the left-hand sides equal to the left-hand sides of the 
equations found in the first order approximation. The 
right-hand sides of the equations are expressed in terms 
of variables of lower order approximations. In carrying 
out this process we obtain that 



U, 



in) 
I 



0(x- n+1 ), t/[ n) = 0(x- n+1 ), n>2. (38) 



This means that the mean velocity starts to behave sin- 
gularly from the third order approximation. 

Taking into account Eqs. (|3"T|) and (|3"5)) we write the 
expansion for the mean velocity in the form 
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(35) 



Since Eq. (|28p has regular singular points, the solutions 
can be obtained in terms of Frobenius series with respect 
to x (for details see, e.g.. Clack et al. [27}, Rudcrman et 
al. 0, Ballai et al. 0). From Eqs. l[3"l" ]) -([3"5" ]) . we see 
that the quantity ^Jj 1 ^ is regular, while all other quantities 

are singular. The quantities vP-\ b£\ and 



behave as In |a;|, while v\J and b L ' behave as x^ 1 , so 
they are the most singular. 

Carrying out calculations on Eqs. ([24]) and (|25|) . and 
utilising Eqs. (|2"g]l. ([HT ]) — ([55 ]> . we find that in the first 
and second order approximations we have 
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d 2 U, 



d 2 U t 
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The functions U^\x), U^'(x), U^'(x) and U\_ LI (x) are 
all continuous at i = 0. Since Eq. (|36[) implies that 
U^\x), U ( f \x), U^(x) and U^\x) are linear functions 



d 2 uf ] 

dx 2 



c?x 2 



(36) 
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of x, we can include U^\x) and U±\x) into [7,| (x) and 

t/| 0) (ir) and take U^\x) = U^\x) = 0, without loss of 
generality. We choose a mobile coordinate system such 
that [7i| (0) = uf \o) = 0. Using this assumption in 
conjunction with Eq. (|3Ci[) it follows that 



lf> = V±x, 



(37) 



where subscript V represents either the subscript '||' or 
'_L' and the functions Vi(x) and V^ n \x) have finite limits 
at \x\ — > 0. The most important property of expansion 
(1391) is that the term of the lowest order approximation 
(proportional to e 1 / 2 ) is very small inside the dissipa- 
tive layer, but becomes large far away from the reso- 
nance. It is also interesting to note that the remaining 
terms tend to zero far from resonance. This result is 
in complete a gree ment with the studies by Ruderman et 
al. (author? ) [26| and Ballai et aL(author?) [37|], where 
they found the mean flow generated by resonant slow 
waves outside the dissipative layer to be proportional to 

,1/2 



IV. INNER SOLUTION 

In this section we determine the inner expansion, which 
is the solution inside the dissipative layer. The thick- 
ness of the dissipative layer is of the order of Zi n hP~ 1//3 , 
where Zi n h is the characteristic scale of inhomogeneity. 
Since we assume that R — GieT 3 ! 2 ) we have l-^Rr 1 ^ — 
^(e^^inh)- As a consequence it is convenient to intro- 

x inside the dissipa- 
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duce a stretching variable £ 
tive layer. 

We can rewrite Eqs. (fT5")) - P3"|) using the stretching 
variable, however, to avoid repetition of previous stud- 
ies and for brevity, we only display the equations we use 
explicitly, which are the perpendicular component of mo- 
mentum and normal component of induction, 
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We use the stretched versions of Eqs. (fl"5"j) - P3"l) to find 
the relationships between variables in each successive or- 
der of approximation (for full details we refer to Clack 
et al. (author?) [27j]). Additionally, Eqs. gDJ) and (gTJ) 
are used to help derive the equations governing the mean 
shear flow inside the dissipative layer. Equations l]24p 
and (|25p for the generated mean flow are transformed to 
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in the first order approximation. This system parallels 
the results of linear theory (see, e.g., Clack et al. [27|). 
For the purposes of this paper, we are only interested in 
the expressions of variables found in terms of vP~\ 
and PW; 
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Following a similar procedure as introduced by Rud- 
erman et aZ.(author?) [26[ we find that variables inside 
the dissipative layer expand in a different way. The large 
variables (v±, b±) have expansion of the form 

g = e^+eg^ + ... ) (44) 

whereas small variables (P, p, p, u, vu, b Xl bu) expand as 

h = eh^ + £3/2^(2) + .... (45) 



The variables vu and vj_ are expanded in series of the 
same form as vu and v±, hence we expand Uu and U± as 



Uu = eUP + 



U± = eU 



(2) 



(46) 



We have chosen the superscripts in Eq. (|4G|) in such a 
way that 



(n) 



rrlll) 



(47) 



On substituting these expansions into the stretched 
versions of Eqs. (fT5"]) -- (f2"3")) we find a linear homogeneous 
system of equations for the variables with superscript l V 



d P w> 

09 



1 rfpw 

d9 



J Aa 



V 



dpo 
dx 



IT [-£ ■ (52) 



In addition the equation that relates and v^P is 



0u^ Ov 



(i) 



99 



■ since. 



(53) 



Here the subscript 'a' means that the equilibrium quanti- 
ties have been evaluated at the resonant position. In the 
second order approximation, after eliminating all vari- 
ables with superscript '1' in favour of and P^ 1 -* using 



Eqs. (|48p (I53p . and satisfying the compatibility condi- 
tion, we derive the equation governing resonant Alfven 
waves inside the dissipative layer(author?) [27j 



(i) 



2jr(i) 



99 



+ V(v + r)) 



v 



V sin a dP^ 



9^ 



POa 



dO 



(54) 



where A a = — (dv 2 A /dx) a cos 2 a, is the gradient of the 
square of the Alfven speed. This governing equation is 
linear despite it being derived by using the full nonlin- 
ear MHD equations, a result explained fully by Clack 
et al. (author?) [23]. Let us summarise their conclu- 
sions needed here for the sake of clarity. If all pertur- 
bations are expanded asymptotically, one can see that 
the second order terms describe magnetoacoustic modes 
and are, therefore, not resonant at the Alfven resonance. 
In addition, these terms act to cancel the small pertur- 
bations produced by the first order modes. This is in 
stark contrast to resonant slow waves, where the second 
order terms act to enhance the perturbations from the 
first order. 



7 



Inserting the expansions (|44 |) -(|46 |) into Eqs. (|42])-([43 
and collecting terms proportional to e, we have 



d 2 t/« 

de 



sinQ! Aw dh \\ 

PoPoa \ x 96 



+ sm a ( t> 




(55) 



magnetic field lines. In order to derive the equivalent 
equation for the generated mean flow perpendicular to 
the magnetic field lines we we note that Eqs. (j40| and 
(|41[) in the second order approximation lead to 



£Rr± B Qa cos a 9b^ 
96 HoPoaV 96 

B 0a ( dB 



V0PQaV Aa 



dx 



sin a dpW 
~p 0a V d6 

1 / dfo 

POa \ dx 



v 9 2 v { l ] 



v 9e 

9vf 



96 



(60) 



rf 2 ^ 
d£ 2 

1 



(1) at; _L 

9Z 



9Z 



9t 



9Z 



II 90 



cos a /,n)9b1>\ B 0a cosa I (1) 9b 



in 



+ _V(^) Utg). (56) 

In deriving Eqs. (J5S]) and we have utilized the fact 
that 



i)=( / I) = (lf) = ( 5 S> =n 



(which follows directly from Eq. (fTT|) ) and Eqs. (|12| and 
(fT5)l . Now we use Eqs. (f4"8")) - ([5"3l to eliminate terms on 
the right-hand sides of Eqs. (f55|) and (|56|) in favour of 

Uj/ and to leave 



d 2 U^ 
d£ 2 



Vsina /^mdP (1) 



d9 



(57) 



d 2 U {2) 



(1) 



96' 



(i) 



POPOa 



9£ 



1 / dpo 



POa \dx J a \ 9£ 

P a V Cfo 



?/' ' /■ > (58) 



Averaging the governing equation (|54[) and substitut- 
ing it into Eq. ((57)) we obtain 



d£ 2 ~ ' '' 



(59) 



which constitutes the equation that governs the gener- 
ated mean flow inside the dissipative layer parallel to the 



,(2) 



Vb 



(2) 



77 sin a dv± f dB$ 

dx 



B 0a cosa V 9£ B 0a 
Substituting Eq. (f5l)l into Eq. (|6T))) we obtain 



^ • (61) 



9v ]_ B 0a cos a 9b ± 
~9f^ ~~ 



POPOaV 9£ 

P a V dx 



y 2 (^ + ?/) 



1 / dpo 

Poa V da; 



1 9 



(gg) 



On substitution of Eqs. (|6"Tj) and (f6"2")> into Eq. 
produce 



(62) 



wc 



d£ 2 



7/ sm a 




v 2 (^ + ?7) 

It follows directly from Eqs. ([53j) and ([54]) that 



(63) 



A £ / (i)~(i) 
— - < u y >v 

V 2 



POaV 



u (Dp(D 



and 



p(i) 



v 



J~(i)9 2 v^ 




9e 



d6 (64) 



(65) 



Using Eq. (J63J) along with Eqs. d64|) — ((65J) we find the 
equation governing the generated mean shear flow inside 
the dissipative layer perpendicular to the magnetic field 
lines is 



d 2 U {2) 




a 3~(l) \ 

x dd 



(66) 



In the next section we will use Eqs. (|59|) and ([66]) along 
with the compatibility condition to derive the jumps in 
the derivatives of the parallel and perpendicular compo- 
nents of the mean velocity across the dissipative layer. 
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V. JUMP CONDITIONS 

To derive the jump in the derivative of the parallel 
component of the mean velocity across the dissipative 
layer we remember from Eq. ([571) that are constants 
and the superscripts '-' and '+' refer to x < and x > 0, 
respectively. This means that 

dul 1] dui 1] 

Vj, + - V,r = lim —± lim 



V{v + rj) 




(67) 



In a similar fashion, we find the jump in the derivative of 
the perpendicular component of the mean velocity across 
the dissipative layer to be 



n - v i 



— lim 



dU 



(i) 



C^oo d£ {-+-oo d£ 



(y + rj) sin a 




d£. (68) 



Outside the dissipative layer we have the following ap- 
proximate equalities: 



M 2 U§\x), U±a^ 2 uf\x). (69) 



II-- -|| 

Let us introduce the new dimensionless variables 



a = 5^x = e^5-% 



1 



V 



where S a = [V (V + rj) / (k\ A a |)] 1/3 , is the width of the 
dissipative layer and k = 2tt/L. Equations (67]) and ((68)) 
can be rewritten in these new variables as 



dU\\ 
dx 



dU± 



dx 



|Aqj cos a 
2ttV 



|A a | sing 
2tt77 



d6 

' X 



d6 



da 



da, 



da. 



(70) 



(71) 



These equations are implicit connection formulae, and as 
such, they must be solved in conjunction with Eq. (|54|) 
(to determine the dimensionless variable q). We note that 
Eqs. ([701) and (JTTJ) imply that at a = 0,tt/2 the jumps 
are zero, since at these values of a there is no Alfven 
resonance present. 

In fact, Eqs. ([7U)) and (fTTj) can be solved explic- 
itly, since the governing equation (|5^|) is linear. It has 
been shown by, e.g. Goossens et aL(author?) [33| and 
Erdelyi(author?) [38J, that the solution of the govern- 
ing equation can be found in terms of the so-called F and 
G functions. In order to find the solutions to Eqs. (170)) 
and (|7ip we need to find in order to obtain q. In 



cartesian coordinates, this was recently accomplished by 
Ruderman(author?) [39] where he found that 



iVP since 
<• = ; . , F(a), 



Po<A|A | 



(72) 



where 



F{a) 



exp (i^>crsgn(A Q ) - <f> 3 /3) d<j>. (73) 



Equation (|72| is the cartesian version of Eq. (66) found 
by Goossens et al. (author?) [331 ]. 

Substituting Eq. (72]) into Eqs. (70]) and (71]) leads to 



dU\\ 



dx 



dU± 



dx 



-A cos a P 2 d9 
Jo 



Asina / P 2 d9 



dF 



da 



dF 



da 



da, 



da, 



(74) 



(75) 



with A = fc 2 sin 2 a/ (27r^pQ a |A Q |). It is shown in the 
Appendix that 



dx 



POO 


dF 


J —OG 


da 


o 


and 


k 2 


sin 2 



da = ir, 



dU ± 
dx 



2^JA Q | 



k 2 sin 3 a 
2^/4 1 A Q | 



P 2 d9, 



P 2 dO. 



(76) 



(77) 



(78) 



Equations ([77]) and ([78]) are the explicit connection for- 
mulae for the jumps in the derivatives of the mean shear 
flow across the dissipative layer. They are explicit be- 
cause we are considering a driven problem, and hence P 
is assumed to be known. 

If we take a = ir/4 we have the following approxima- 
tion 



IdUf 




'dU±_~ 


dx 




dx 



! V2y 

7. ' 
1 1 n n 



(79) 



and these values can be seen as jumps in vorticity. Here 
we have used the obvious estimates |A a | = 6{y 2 /Zi n h), 
P = ff{€V 2 khnb) and V = /ff(e 3 / 2 Vl inh ). In order to find 
the profiles of the components of the generated mean 
shear flow we need to impose boundary conditions far 
away from the dissipative layer. For example, if there 
are rigid walls at x = ±o where the condition of adhesion 
has to be satisfied, then the components of the generated 
mean flow take the simple form 



Un = 



dx 



dx 



2 ' 
x + a 



x > 0, 



, x < 0, 



(80) 
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U ± = 



dx 
dU 



dx 



2 

x + a 



, x > 0, 



(81) 



x < 0. 



With the estimate given by Eq. ([79")) and the simple 
mean flow profiles in Eqs. (|80[) and (|STj) we can calculate 
the expected mean flow generated outside the dissipative 
layer in the solar upper chromosphere and solar corona. 
For example, if the incoming wave has a dimensionless 
amplitude of e = i^(1CP 4 ), then the predicted mean shear 
flow is of the order of lOkms -1 in both the upper chro- 
mosphere and corona. Here we have assumed that the 
characteristic scale of inhomogeneity (/j n h) is 10 2 m in the 
upper chromosphere and 10 3 m in the corona. This gen- 
erated flow can be superimposed on existing flow, so it 
is difficult to observe. These results should be used with 
caution, when applied to the solar atmosphere. We have 
assumed that outside the dissipative layer the plasma is 
homogeneous and infinite, and clearly the solar atmo- 
sphere is neither, so, since all flows are entirely governed 
by the boundary conditions, the flow generated may be 
less than predicted or may even be stopped entirely. 



VI. CONCLUSIONS 

In the present paper we have completed the nonlin- 
ear theory of resonant Alfven waves in dissipative layers 
in a one-dimensional (1-D) planar geometry. Clack et 
al. (author?) [27j showed that even though nonlinear- 
ity and dispersion are considered the equation governing 
resonant Alfven waves in the dissipative layer is always 
linear (provided e <C Rr 1 ^). However, they neglected 
the second manifestation of nonlincarity at resonance; 
the generation of mean shear flows outside the dissipa- 
tive layer. This flow is produced by the nonlinear inter- 
action of harmonics inside the dissipative layer and may 
still exist even though the governing equation inside the 
dissipative layer is linear. 

We have shown that outside the dissipative layer a 
mean flow is generated parallel to the magnetic surfaces. 
The flow has an amplitude proportional to e 1 / 2 , and de- 
pends linearly on x. The derivatives of the velocity of the 
generated mean flow have a nonzero jump across the dis- 
sipative layer determined by Eqs. (|70|) and (fTTj) . When 
a = it/ 4 the magnitude of the jumps can be estimated 
as e 1 / 2 ^"^, where V is the phase speed of the incoming 
wave and Z; n h is the characteristic scale of inhomogeneity. 
For typical conditions in the solar upper chromosphere 
and corona the magnitude of the jumps in the deriva- 
tives of mean shear flow velocity would be of the order 
of 10s _1 . From this, and the simple flow profiles given 
by Eqs. ([80)1 and (|8Tj) , we predict a mean flow outside 
the dissipative layer (generated by resonant absorption) 
with an amplitude of the order of lOkms -1 in the solar 
upper chromosphere and corona. 



The magnitude of the jumps in the derivatives of the 
mean flow were found to depend on the phase speed (V), 
the characteristic scale of inhomogeneity (^ n h) and the 
dimensionless amplitude of oscillation (e) far away from 
the dissipative layer. Even though simple estimations 
allowed us to approximate the magnitude of the shear 
flows, the results in the present paper should be used 
with caution. The present analysis has been carried out 
for magnetic configurations that are homogeneous and 
infinite in the direction of wave propagation outside the 
dissipative layer. This situation can only take place in 
laboratory devices (such as tokamaks). In the solar atmo- 
sphere, magnetic configurations are bounded and/or in- 
homogeneous in the direction of wave propagation. These 
additional boundary conditions may reduce (or even pre- 
vent) the generation of mean shear flows by resonant ab- 
sorption. 

It is worth mentioning that when an Alfven resonance 
is present there may exist a slow resonance coupled to 
it. If this happens, both resonances will generated their 
own mean shear flows which could become superimposed. 
This may enhance the flow through constructive superpo- 
sition, but equally, they may interact destructively leav- 
ing a reduced overall combined generated mean shear 
flow. In other words, it is possible to see a more effi- 
cient coupled dissipative layer create more heating, but 
actually generate a smaller amplitude mean shear flow 
compared with two separate dissipative layers. 

If the boundary conditions allow the generation of a 
mean shear flows by resonant absorption in the solar at- 
mosphere, we would expect to see these flows superim- 
posed on existing large scale motions. It is speculated 
that these mean shear flows (if set up) superimposed on 
existing bulk motion of the plasma would produce tur- 
bulent motions and create a Kelvin Helmholtz instability 
(KHI) along the dissipative layer. A KHI instability may 
already exist due to the shear velocity field of Alfven 
waves as was shown by Terradas et al. (author?) [40] |. 
The generation of mean shear flow would survive the in- 
stability due to the velocity shear and, hence, will still 
contribute to the bulk motion causing larger shear, and 
perhaps creating a secondary instability. Either way, the 
KHI would distort the dissipative layer, locally enhanc- 
ing the dissipation in (due to enhanced transport coeffi- 
cients), which in turn would produce greater amplitude 
mean shear flows. This process of feedback would con- 
tinue until either; (i) the wave is dissipated, or (ii) the 
turbulence destroys the dissipative layer. Either way, the 
turbulence would allow greater absorption and achieve 
heating over a greater area (here the development of tur- 
bulent eddies can be seen as a means of heat transport). 

The resonant absorption itself may not produce the 
heating we see in the solar corona, but we believe it is 
likely that these dissipative layers are created and KHI 
are formed. The KHI might be able to trigger small scale 
reconnection events (nanoflaring) if the mean shear flows 
produced by resonant absorption become strong enough. 
However, further study and observations are required be- 
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fore we can make more definitive investigations. 
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are calculating the absolute value of the derivative, and 
so we use complex conjugates. If we change the order the 
integration we obtain 



^e"^/ 3 d<f> 



poo 


dF 


2 p 




da = 


' — oo 


da 


Jo 






roc 




X 








Jo 



The third integral in Eq. (|A2[) is the definition of the 
delta function, S(X ~ (f>), so we can write Eq. (|A2|) as 



APPENDIX: THE CALCULATION OF THE 
INTEGRAL OF THE F FUNCTION 

In this Appendix, we will derive Eq. (|76p . We know 
the form of F(a) from Eq. (|T5|) , so we can rewrite Eq. 
||75) as 



dF 



da 



da 



da 



o i0a-0 3 /3 



II 



I Xe~ lXa - x /3 dX ) , (Al) 



where the tilde denotes the inclusion of the sgn(A a ), from 
this point onwards we drop the tilde notation. Here we 



f 

J — i 



dF 



da 



da 



0e^ 3/3 dej) 



x / Ae~ A3/3 • 27n5(A - <f>)d\. (A3) 



Since the integral of the delta function is always unity, 
we arrive at 



poo 


dF 






da 


da = 2ir [ 


' — oo 


Jo 



-7TC 



-2<£ 3 /3 



(A4) 
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